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Abstract 

The Laplace transforms of the transition probability density and distribution functions for the 
Ornstein-Uhlenbeck process contain the product of two parabolic cylinder functions, namely 
Dy (x) Dy (y) and Dy {x) -D„_i (y), respectively. The inverse transforms of these products have 
as yet not been documented. However, the transition density and distribution functions can 
be obtained by alternatively applying Doob’s transform to the Kolmogorov equation and cast¬ 
ing the problem in terms of Brownian motion. Linking the resulting transition density and 
distribution functions to their Laplace transforms then specifies the inverse transforms to the 
aforementioned products of parabolic cylinder functions. These two results, the recurrence rela¬ 
tion of the parabolic cylinder function and the properties of the Laplace transform then enable 
the calculation of inverse transforms also for countless other combinations in the orders of the 
parabolic cylinder functions such as Dy (x) n „_2 (y), (x) n„_i (y) and Dy (x) n „_3 (y). 
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1 Introduction 


The parabolic cylinder function is intensively used in, for instance, chemical physics (Kalmykov, 
Coffey and Waldron m), lattice field theory (deLyra, Foong and Gallivant [2]) and astrophysics 
(Zaqarashvili and Murawski ^). This paper obtains the inverse Laplace transforms for products 
of two parabolic cylinder functions by exploiting the link between the parabolic cylinder function 
and the transition probabilities of the Ornstein-Uhlenbeck process that is used in, for example, 
astrophysics (Chandrasekhar [3]), neurophysiology (Gluss [S]) and hnance (Vasicek [B]). 

The Laplace transforms of the transition probability density function - in short, the transition 
density - and the transition probability distribution function - in short, the transition distribution 
- of the Ornstein-Uhlenbeck process contain products of two parabolic cylinder functions for which 
the Laplace-transform parameter appears in the order of both functions. Their inverse transforms 
have as yet not been documented. However, the transition density and the transition distribution 
of the Ornstein-Uhlenbeck process can be obtained by applying Doob’s transform to the Kol¬ 
mogorov equation (Cox and Miller |5], Risken [S]) in view of casting the latter equation in terms 
of Brownian motion for which solutions are well-documented. Combining the resulting transition 
density and transition distribution with their Laplace transforms then yields the desired inverse 
transforms for Laplace transforms that contain, in stylised form, [x) (y) (transition density) 

and Dy {x) (y) (transition distribution). The latter two results, the recurrence relation for the 

parabolic cylinder function and the properties of the Laplace transform then also enable the calcu¬ 
lation of inverse transforms for countless other combinations in the orders of the parabolic cylinder 
functions, such as (x) (y), (x) (y) and (x) (y). These results are new to 

the literature as the inverse Laplace transforms for products of two parabolic cylinder functions in 
Erdelyi et al. m and Prudnikov, Brychkov and Marichev m do not allow the Laplace-transform 
parameter to emerge in the order of the parabolic cylinder functions. 

The inverse Laplace transforms that are derived in this paper can be specialised in interesting 
directions. For instance, employing x = y can yield expressions for the inverse transform of squared 
parabolic cylinder functions. Setting either x = 0 or y = 0 returns expressions for single parabolic 
cylinder functions and the limiting case of x = y = 0 produces inverse transforms for ratios of 
gamma functions. 

The remainder of the paper is organised as follows. Section 2 lists the properties of the Laplace 
transform and the parabolic cylinder function that will be used in the paper. Section 3 derives 
the Laplace transforms and the transition density and transition distribution for the Ornstein- 
Uhlenbeck process. Section 4 illustrates how the resulting two inverse Laplace transforms can be 
extended to other orders. Section 5 presents four specialisations for single parabolic cylinder func¬ 
tions and ratios of gamma functions and discusses their relation to expressions that are documented 
in the literature. 
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2 Properties of the Laplace transform and the parabolic cylinder 
function 


Extensive detail on the Laplace transform and the parabolic cylinder function can be found in, for 
instance, Debnath and Bhatta [12] and Erdelyi et al. m, respectively. This section therefore only 
lists the properties that are employed in deriving the below inverse Laplace transforms. 

2.1 The Laplace transform and its properties 

The Laplace transform of the original function / (t) is 

H-CO 

L {/ (t)} = 7 (s) = j exp {-St) f (t) dt, 

0 

where Res > 0. The inverse Laplace transform is denoted by 

L-m{s)} = f{t). 

The linearity property of the Laplace transform gives 

L {af (t) + bg (f)} = of (s) + bg (s). (2.1) 

Differentiation of the original function / (t) yields 

L{/'(f)} = s7(s)-/(0) (2.2) 

and linear transformation in the image function gives 

L{exp(af)/(f)} =7(s-a). (2.3) 


2.2 The parabolic cylinder fnnction and its properties 


The parabolic cylinder function of order v with argument z is denoted by {z). Equation (8.2.9) 
in Erdelyi et al. m specifies the parabolic cylinder function for n = 0 as 

(2.4) 


Do (z) = exp 1^- 
Equations (19.3.1) and (19.3.5) in Abramowitz and Stegun 


give 


1 — V 


(2.5) 


where T (v) denotes the gamma function. Simplifications in the below results often employ the 
recurrence and the duplication formulas for the gamma function 

T (n + 1) = uT (v ), 


3 


(2.6a) 

(2.6b) 





see Equations (6.1.15) and (6.1.18) in Abramowitz and Stegun [14| . 

The derivative of (z) to its argument is given by 

D'v (z) = (z) - Dy+i (z), (2.7a) 

n; (-z) = ^zD, (-z) + (-z), (2.7b) 

as can be seen from Equation (8.2.16) in Erdelyi et al. [T3j. Equation (8.2.10) of Erdelyi et al. [T3j 
specifies the Wronskian as 

D, (z) n; (-Z) - D, i-z) n; (z) = (2.8) 

T{-v) 

The recurrence relation for the parabolic cylinder function is given in Equation (8.2.14) of Erdelyi 
et al. [T^ 

Dy+i (z) - zDy (z) + vDy_i (z) = 0. (2.9) 

The below identity directly follows from Equations ()2.7al) - (j2.9p 


{z) (-z) + Dy {-z) (z) = 




—uT {—v )' 

Also the following two integral expressions can be obtained from the above properties 


( 2 . 10 ) 


I^v—1 (■2^) ? 

(2.11a) 

Dy_l (-Z) . 

(2.11b) 


3 The transition density and transition distribution of the Ornstein- 
Uhlenbeck process and their (inverse) Laplace transform 

Let W = {1L7, t ^ 0} be an Ornstein-Uhlenbeck process with initial value Wq = wq. Its dynamics 
are given by 

dWt = {a - ^Wt) dt + adZt, (3.1) 

with jd > 0 and t > 0. The positive value of /3 implies that this stochastic process reverts to its 
long-term mean that is at The instantaneous variance is given by cr^, with ct > 0, and dZt is 
the increment of a Wiener process. 

The transition distribution is defined as P {w, t\wo) = Pr {W (t) ^ r/;|iy (0) = rco} and the tran- 

d 

sition density then is p{w,t\wo) = -^P {w,t\wo). The transition density specifies the probability 
of attaining w at time t given that the process initially is at the source point wq. The transition 
density satisfies the Kolmogorov forward and backward equations as discussed in, for instance, Cox 
and Miller [5] and Risken . The Kolmogorov backward equation for the process in Equation (|3.1I) 
is 


1 2 d‘^Piw,t\wo) dp{w,t\wo) dp{w,t\wo) 

T — ftS — + ~ au,o =— m —■ 
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The initial condition is given by 


p {w, t|tco) = S {w — wq) for t = 0, 

where (5 (•) is the Dirac delta fnnction that gnarantees that all initial probability mass is located 
on the initial state of the process. 

The Laplace transform of the transition density, p{'w,s\wo), is 


+OD 

p{w,s\wo) = / exp {—st)p{w,t\wo) dt, 


with Res > 0. The Kolmogorov backward equation then can be rewritten as 
1 2 d‘^Piu!,s\wo) , ^ ^ ^dp{w,s\wo) I , 

2^ --^ “ PWo) -- sp{w, sl^o) = -d{w- rco) • (3.2) 

The substitutions zq = wq—^ and go = and the transformation v {q, s|go) = exp P (?) ^ko) 

simplify the homogenous part of Equation (j3.2p into 


d‘^v {q,s\qo) 

dql 


+ {^-^-^\viQ,s\qo)=0 


The parabolic cylinder functions (^o) and (—Qo) are two linearly independent solutions 

to the latter equation (see Equation (8.2.5) in Erdelyi et al. [IS]). The resulting Green’s function, 
G{w,wo), is 


G{w,wo) = G{w,wo)\_ 


OO^WQ^W 


+ G{w,wo)\ 


W^W0^-\-OO ’ 


. , „, M . f i/dwo — a)^ 

with G{w,wo)\_^^^^^^ = Aexpl —— 


D 


-s/p 


G{w,wo)\ 


W^W0^-\-OO 


= B exp 


(Pwq - py 


D 


-s/p 


( V2 {(3wo - a) 

\ CTV^ 

f V2{j3wo - a;)\ 




cr 


VP 


for which the constants A and B are to be obtained from the conditions 


G(u^,«^o)L^^o^+oo - = 0 for wq = w, 

dG(in,rco)Ls:^o^+oo dG 2 

- ;--;-- =-for Wq = W. 

dwQ dwo 

The derivatives and the Wronskian in Equations (|2.7ap ~ (j2.8p allow to express the Laplace transform 
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of the transition density as 


p{w,s\wo) 


( 


< 


r {sj^) I {f3wo - af {/3w - a)^ 
I 2a2/3 


xD 

r(V/3) 

xD 


-s/p 

exp 

-s/p 


fV2{l3w-a)\ ( V2{f3wo-a)\ 

\ i 

/ (/3wo — a)'^ {/Pw — a)'^\ 

I 2^ 2cj2/? j 

( V2{Pw-a)\ fV2{Pwo-a)\ 


for — oo ^ Wo ^ w, 


for w ^ wq ^ + 00 . 


(3.3) 

However, the literature did not obtain the original function p(w,t|wo) by inverting the Laplace 
transform (|3.3p . Instead, Doob’s transform, see Doob [7], offers a simpler route by simplifying the 
Kolmogorov forward equation into the diffusion equation for which solutions are widely available. 
The resulting transition density p{w,t\wo) is documented in, amongst others, Gluss [5] and Cox 
and Miller i) as 


p{w,t\wo) 


7rcj2 (1 — exp (—2/3f)) 



((/3w 


— a) — (/3wo — a) exp (—/3f))^ 
/3 (t 2 (1 — exp {—2/3t)) 


(3.4) 


The transition density (|3.4p thus specifies the inverse Laplace transform of Equation (|3.3p . 

The transition distribution, P {wi,t\wo), specifies the cumulative probability mass that is situ¬ 
ated below wi- Its Laplace transform is 


P(wi,s|wo) 


+OD 

J exp (—st) P {wi,t\wo) dt, 
0 


where Res > 0. The Laplace transform P(wi,s|wo) can be obtained from the Laplace transform 
of the transition density by employing the definition of the transition distribution as follows 

+ 0O f Ull 

P(wi,s|wo) = / exp(—st)< J p (w, f|wo) dw > dt 


0 

( +00 


P(wi,s|wo) = 


exp {—st) p {w, f|Wo) dt > dw 


—oo V 0 

Wl 


p 


(wi,s|wo)=y p(w,s|wo)dw. 


The expression for the Laplace transform P (wi, s|wo) will differ in function of the location of wi 
with respect to wq. However, the precise position of wi does not affect the below results such that 
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we continue here with the restriction wi ^ wq. The Laplace transform (j3.3p gives 


P{wi,s\wo) 


r(s//3) 


exp 


/ {Pwq - af\ 
\ 2ct2/3 ) 


D-sip 


V2{Pwo - 

a^/^ J 


X 



(/3w — a)^\ 

2 ^ 2/3 J 


^-s/P 


( V2{/3w - «) \ 

\ cjv^ j 


dw 


+ 



V2{/3wo - a) \ 
crv^ j 


y/2 {/3w — a) \ 

crVP J 


dw. 


The integrals (j2.11ap and ()2.11bh and the identity ()2.inp then allow to specify P {wi, slrco) as 

{(3wo-af {I5wi-af\ 


-P! 1^1 r(s//3) 

P(u;i,s|u;o) =- exp 

s /3v27r 


X D 


-s/d 


2a^P 

V2{/3wo — a 




D 


-i-s/d 


2a‘^P J 

( \/2(/3wi - a) 

\ CJv^ 


(3.5) 


for wi ^ wq. 


Equation (j3.5p has a clear probabilistic interpretation as the inverse transform of the term — is 
1. The latter term thus corresponds with the cumulative probability of the Ornstein-Uhlenbeck 
process over its entire domain, i.e. (—oo, Too). The second term then must be seen as representing 
the Laplace transform of the loss of probability that arises when wi takes on a value below +oo. As 
required, lifting tci to +oo reduces the second term to zero as can be seen from Equations (19.3.1) 
and (19.8.1) in Abramowitz and Stegun |14| . 

The original function P {wi, s|t(;o) can directly be obtained by integrating the transition density 

dial) 


P (tcijflrco) 


-erfc — 


{{Pwi - a) - {Pwq - a) exp j-fdt)) 
A//3cr2(l -exp(-2/3t)) 


(3.6) 


where erfc( 2 ;) denotes the complementary error function, see Abramowitz and Stegun |14j . 

The inverse transform of the Laplace transform (13.511 thus is given by the transition distribution 

(IMI). 


4 Inverse transforms of Laplace transforms that contain products 
of two parabolic cylinder functions 

This section illustrates how the recurrence relation of the parabolic cylinder function and the above 
properties of the Laplace transform can be used to obtain further results out of the two central 
inverse Laplace transforms. However, the below 10 inverse transforms, that are also collected 
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in Table 1, thus must be seen as a subset of the infinitely many inverse transforms that can be 
obtained. 


Table 1: around here. 

We start with the (inverse) Laplace transform of the transition density in Equations (j3.3p and (|3.4p 
for —oo ^ Wo ^ w. Notation is simplihed by using x = — and y = The 

linearity property (12.11) allows to re-express the Laplace transform in Equation (13.311 as /(s) = 
T (s//3) D_gip (x) D_gijj (y) for which the previous condition — oo ^ wq ^ w translates into x + y 
0. Note that the latter condition does not impose any additional (linear) relation between the two 
arguments. Linearly transforming the image function via Equation (|2.3I) for s-|-c with c ^ 0 creates 
a more versatile specification for the order of the parabolic cylinder functions 

L-i {r (^) ( x ) Il-(.+e)//3 (y)} = 

/3exp(-ct) ^(y^-x‘^\ _( (y-hxexp(-/3t))^^ 

- ^l-exp(-2/3t)"''^V 4 2(l-exp(-2/3t))j 

[Re s > 0, /? > 0, c ^ 0, X -|- y ^ 0] . (4.1) 

It is to be noted that Equation (4.1), with the same condition, also arises when evaluating the 
Laplace transform in Equation (|3.3I) for the domain w ^ wq ^ -(-oo. 

Applying the same steps to the (inverse) Laplace transform for the transition distribution in 
Equations (|3.5p and (13.61) gives 

L-i {r (^) (x) ii_i_(.+e)//3 (y)} 

^ P^/^^exp{-ct) / y^ - x^ \ / yxexp (-/3t) \ 

a/2 4 J y y/2 (1 - exp (-2/3t)) j 

[Re s > 0, /3 > 0, c ^ 0, X -|- y ^ 0] . (4.2) 

Fixing the order of the hrst parabolic cylinder function in the product at — {s + c) /f3 and using 
the recurrence relation (US]) on the second parabolic cylinder function with v = —1 — {s + c) / P 
gives 

£±|±£r (^) (x) I?- 2 -(«+c )//3 (y) = (4.3) 

~ ( /3 ) ^—(s+c)/y (x) L1_i_(s-)-c)//3 (y) + T ^ ^ ^ D—(^s+c)/i3 ( x ) H_(s-(-c)//3 (y) ■ 

Plugging the inverse Laplace transforms (SH) and (14.21) into relation p4.3ll gives 

L-i {£±|±2r (^) l?_(.+e)/y (^) ^- 2 -(.+c)/y (y)} 

/3exp(-ct) f (y + xexp(-/^t)y'] 

4 J 1^1 _ exp (-2/31)"''^ 2(l-exp(-2/3t))j 

y/3y/7fexp(-ct) j yx exp (-/3t) ^ | 

V2 \a/ 2(1 -exp(-2/3t)) j j 

[Re s > 0, /3 > 0, c ^ 0, X -|- y ^ 0] . (4.4) 
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Increasing n to — (s + c) //3 in the recurrence relation (12.91) yields 


^ ^ ^ ^—(s+c)/fi (®) ^i—(s+c)/f} iy) 


(4.5) 


yT (^) Z?_(,+e)/^ (x) {y) + (^) (x) I)-l-(.+c)//3 iv) ■ 

The inverse transform of sT ( D_(^s_^_cy^ (x) [y) is obtained by applying the differ¬ 


entiation property (j2.2l) to the inverse transform (|4.2p . The term / (0) is for x -|- y = 0 and 0 
otherwise, for which we employ the notation . Indeed, the complementary error function 


x+y=0 


for X -)- y ^ 0 and t = 0 gives erfc(-l-oo) = 0, whereas it yields erfc(O) = 1 for x -|- y = 0. The inverse 
transform then can be written as 


L 


-1 


{sr (if) (!■) (y) - 


= exp 


y‘^-x^\ I exp (- (/3 c) t) 
[ (l-exp(-2/3t))i 


(x-hyexp(-^t))exp - 


(y-Hxexp(-/3t))^ 
2 (I — exp (—2/3t)) 


c/3v^exp (-ct) „ ( y-I-X exp (-/It) 

-;=-erfc — , 


V2 


\\/2(I -exp(-2/3t)) 


[Re s > 0, /? > 0, c ^ 0, X -|- y ^ 0] . 
Plugging the inverse transforms (|4.1D . (j4.2p and 


(4.6) 


into relation (|4.5p yields 


L 


-1 


( y ) ^-is+c )/0 (x) T>i_(s+c)//3 (y) 
j3 exp (—ct) 


(y -I- xexp {-fit)) exp 


(I-exp(-2/3t))2 
[Re s > 0, /I > 0, c ^ 0, X -|- y ^ 0] . 


x+y=Q 

y^ — x^ 


exp 


(y-Hxexp(-^t))^ 
2(1 — exp (—2/3t)) 


(4.7) 


Fixing the order in the second parabolic cylinder function at —1 — {s + c) /fi and using v = 
— (s -|- c) /fi within the recurrence relation (12.9D for the first parabolic cylinder function gives 

r (^) T»1 -(s+c)//3 (x) T»_i_(s+c)// 3 (y) = (4.8) 

xT t>_(s+c)//3 (x) T)_i_(^+c)// 3 (y) + t>_i_(^+c)//3 (x) t>_i_(s+c)// 3 (y) • 

The inverse transform of ^^F -^-i-(s-i-c)//3 (®) ^-i-{s+c)/p iv) is obtained via applying the 

linear transformation p2.3p for f {s + ji) to Equation (14.111 and using the recurrence formula for the 
gamma function in Equation (I2.6bh 

I)-i-(.+c)/y (x) I)-i-(.+c)//3 (y)} 

_/3exp(-(/3-Fc)t) ^(y‘^-x^\ ( (y-h xexp (-/It))^^ 

“ Vl - exp (-2/3t) 4 J [ 2 ( 1 - exp (-2^t ))) 

[Re s > 0, /I > 0, c ^ 0, X -|- y ^ 0] . (4-9) 
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Connecting the inverse transforms (|4.2D and (14.911 to the recurrence relation (14.81) yields 


L-l |r ^l-(s+c)//3 (a^) I)-l-(.+c)//3 I 

= ex f / /3exp(-(/3 + c) t) 

^ V 4 y I - exp (-2/3t) 



(y + xexp(-/3t)f '\ 
2 (1 — exp (—2/3t)) y 


x/30rexp (-ct) „ 

H-=-eric 

V2 


/ 7/+ xexp (-/3t) 

V\/2(l -exp(-2^t))^ 


[Re s > 0, /3 > 0, c ^ 0, X + 7/ ^ 0] . 


(4.10) 


Setting the order in the second parabolic cylinder function at t; = —2 — (s + c) //3 and keeping the 
order of the first one at — (s + c) //3 gives 

£±|±£r (x) 2^_3-(.+c)//3 (?/) = (4.11) 

~ 2/4^ ^ ^—{s+c)/p (a^) ^—2—(s+c)/p iu) 4" 4^ (a^) 44_i_(s+c)//3 (l/) ■ 

The inverse transform for the hrst term on the right-hand side of relation (|4.11l) can be obtained 
from earlier results in two steps. First, applying the linear transformation in Equation (|2.,'ip for 
f {s + f3) and the recurrence formula for the gamma function (|2.6ap to the inverse transform (14.101) 
gives 


L 




D 


-2-{s+c)/P 


= exp 


7/^ — ' 


/3 exp (— (2/3 -I- c) t) 
V^l -exp(-2/3t) 


■ exp 


{y + X exp I3t)) 
2(1 — exp (—2/3t)) 


^ x/jy/Trexp {-{13 + c)t) / y -F xexp {-fit) 

y V^2 (1 - exp (-2/3t)) 
[Re s > 0, /3 > 0, c ^ 0, X -|- 7/ ^ 0] . 


(4.12) 


Second, subtracting the latter result from the inverse transform (14.411 gives 

L-i |r I4-(.+.)//3 (x) 44-2-(.+c)//3 (y)} 

(y‘^ - x‘^\ \ a r-i - 1 I (y+ a;exp(-/3t))^^ 

= exp J |/3V4-exp(-2/3t)exp (-cf) exp ^ (i _ ^xp (-2/3t )) ) 

(xexp {-(3t) + y) /j^/Trexp (-ct) ^ y + xexpj-fit) ^ | 

V\/2(l -exp(-2/3t)) y J 

[Re s > 0, /3 > 0, c ^ 0, X -|- 7/ ^ 0] . (4-13) 


Finally, combining the inverse transforms (|4.2p and (14.131) with the recurrence relation (14.1111 yields 
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the desired result 


L-" |£±|±£r (^) z)_(,+,)/;3 {x) i?-3-(.+c)//3 (y)} 

(y‘^-x‘^\\ ^ ^ ( {y + xex.p{-pt)f 

= exp |-!/;3Vl-exp(-2,3«)exp(-d)exp 

(?/(xexp(-/3t)+ y) + l)/30rexp(-ct) [ y + xexp(-/3t) \1 

^-7!-“* i75?T=SFF5lj I 

[Re s > 0, /? > 0, c ^ 0, X + ?/ ^ 0] . 


(4.14) 


5 Some specialised results 

The above inverse Laplace transforms can be specialised by equating the arguments (in absolute 
value) and/or by setting one or both of them at 0. This section briefly specifies four such speciali¬ 
sations. 

First, an inverse Laplace transform for squared parabolic cylinder functions emerges from the 
result in Equation (14.11) by setting y = x 

/3exp(—ct) / (1-b exp (—/3t))^\ 

y^l — exp {—2fjt) ^ 2 (“2/3t)) J 

[Re s > 0, /3 > 0, c ^ 0, X ^ 0] . (5.1) 

Second, plugging y = —x into the inverse transform (14.71) gives 

L-I {r (^) 7?_(,+e)//3 {x) (-X) - 

—/3x (1 — exp [—/3t]) exp (—ct) / x^ (1 — exp (—/3t))^\ 

(l-exp(-2/3t))5 2(l-exp(-2/3t)) J 

[Res > 0,/3 > 0,c ^ 0] . (5.2) 

Note that the additional factor ^ is to be included irrespective of the value of the argument as 
the condition x + y = Q for this factor in Equation (14.7p always holds. Also, the inverse transform 
(15.21) now applies to all real values for the argument. As required, both sides in Equation (|5.2I) 
vanish when additionally restricting x to 0 as can be seen from property (12.5p and the duplication 
formula (I2.6bl) . 

Third, the results in the previous section generate a large body of expressions for single parabolic 
cylinder functions. Eor instance, using y = 0 in Equation (|4.ip gives 

^-1 |2(.+c)/(2/3)p I?-(s+c)//3 (X)} 

2(3 exp {—ct) / x^ (1-|-exp (—2/3t)) \ 

^yY^^expJ^^2^ 4 (1 — exp (—2/3t)) y 

[Re s > 0, /3 > 0, c ^ 0, X ^ 0] . (5.3) 
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The result in Equation (j5..SD complements extant inverse transforms for single parabolic cylin¬ 
der functions that are documented in Erdelyi et al. m and Prudnikov, Brychkov and Marichev 
m- For instance, the Laplace transform in Equation (5.18.9) in Erdelyi et al. |10] is / (s) = 
(s-|-u) Zl_ 2 s (x) and its original function is given by f (t) = exp (exp (t) — 

exp ^ D 2 v ^ Laplace transform in (15.811 thus offers a more flex¬ 

ible specification for the order in the parabolic cylinder function. The Laplace transform (5.18.9), 
on the contrary, does not impose the tight connection between the argument in the gamma function 
and the order in the parabolic cylinder function that is present in Equation (15.3p . Note that both 
specifications are identical when employing c = 0 and /? = ^ in (15.311 and specialising (5.18.9) for 
u = 0 and using the property (|2.4I1 . 

Finally, setting both arguments to 0 yields inverse Laplace transforms for ratios of gamma 
functions. For instance, plugging x = y = 0 into Equation (14.1211 gives 



[Re s > 0, /? > 0, c ^ 0] . 


2/3 exp (— (2/3 -|- c) t) 
V^7r(l-exp(-2/3t)) 


(5.4) 


Erdelyi et al. m and Prudnikov, Brychkov and Marichev na list a large number of inverse 
transforms for ratios of gamma functions. The inverse transforms that can be obtained from the 
results in this paper are complementary as the expressions in the literature multiply the Laplace- 
transform parameter either by 1 or 2 (see for instance Sections 3.1.1. and 3.1.2 in Prudnikov, 
Brychkov and Marichev m), but at the same time allow for a more flexible connection between 
the arguments in the various gamma functions. 
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* Re s > 0, /3 > 0, c ^ 0, X + j/ ^ 0. 
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Table 1: continued* 


+ 00 


7 (s) = J exp (-st) / (t) dt 

0 

i"M/w} = /w 


7. 


^-2-(s+c)//3 (y) exp I 


^ ( y + xe^vi-pt) 

i V2(l-exp(-2/3t))^ 


y"^ - x^\ ] P exp (- (2/3 + c) t) 


X exp 


a/I - exp (-2/3t) 

(// + a:exp(-/3t))^ \ 

2 (1 — exp (—2/3t)) y 


+ 


cc/jy/TT exp (- (/? + c) t) ^ y + a:exp (-/3t) 

1^/2 (1 - exp (-2/3t)) 


s+ff+c -p 
o -L 


s + c 


£>_(s+c)// 3 (a;) D_ 2 -(^s+c )/0 (y) exp 


2 2 ' 
y — X 


X exp — 


j /3exp(—ct) 

I a/ 1 - exp (-2/3t) 

(// + xexp(-/3t))^ \ 

2 (1 — exp (—2/3t)) I 


///3A/7rexp (-ct) ^ y + xexp (-/3<) 

a/ 2 yA/2(l-exp(-2/3t))^ 


9. ^±2^r f ] £)_(s+c)// 3 (a;) i:>_ 3 -(s+c )//3 (y) exp ( ^ ^ ) <j -y/3\/l - exp(-2/3t) 


+ 


10 . 


^ ^ ^ ^-l-{s+c)/l3 (x) -D-1 -(s+c)//3 (y) 


X exD (-ct) exp f _ (v+^^^rLM)!] 
xexp( ctjexp^^ 2 (1 - exp (-2/3t)) j 

(y (xexp(-pt) +y) + 1)/SA/irexp (-ct) 

7i 

X erfc (J±£2Eh«_) I 

I A/2(l-exp(-2/3t))y J 


/3exp(-(/3 + c)t) (y'^-x" 
• exp ' 


a/1 - exp (-2/3t) 


X exp — 


{y + xexp{-pt)f \ 
2 (1 — exp (—2/3t)) j 


* Re s > 0, /3 > 0, c ^ 0, X + y ^ 0. 
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